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1. Introduction 

The Heavy Nuclei Experiment on HEAO-3 (Binns et al., 1981) uses a combination 
of ion chambers and a Cerenkov counter. During analysis, each particle is assigned two 
parameters, Zc and Z h proportional to the square roots of the Cerenkov and mean ioniza- 
tion signals respectively. Because the ionization signal is double valued, a unique assign- 
ment of particle charge, Z, is not possible in general. Our previous work (Binns et al., 
1983, 1985, and Stone et al., 1983) has been limited to particles of either high rigidity or 
low energy, for which a unique charge assignment was possible, although those subsets 
contain less than 50% of the total number of particles observed. In this paper we discuss 
the use of the maximum likelihood technique to determine abundances for the complete 
data set from —1.5 to —80 GeV/amu. 

Figure 1 shows the possible values of Zc and Z x for elements near iron, and indicates 
the substantial overlap between adjacent elements, even before smearing by the resolution 
function. In Figure 2, the curves of Figure 1 have been transformed using the variable 
Zc/Zi instead of Z\. This transformation simplifies the following data analysis. 

2. Analysis 

Particles were selected from the full exposure, 580 days, and were required to have a 
good Cerenkov signal, at least one good ion chamber and a reliable trajectory. These par- 
ticles were assigned an initial charge estimate, Z^, and 1/40 of those with Zes,>19.5 were 
saved, together with all the remaining particles with Z est >30. The selected particles were 
binned in a two dimensional histogram, with one axis being the logarithm of Zc and the 
other being Zc/Zi. Figure 3 shows a contour plot of the region of this histogram near 
iron. 
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Figure 1. Curves of Z\ versus Zcfor the elements near iron. 
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Figure 2. The curves of Figure 1, displayed in (Zi/Zc, Zq) space. 
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We have used the iron distribution as a reference distribution for the other elements, 
by scaling it according to the following criteria: 

1) All elements are assumed to have the same energy spectrum, although it is known 
that the sub-iron secondary elements have steeper spectra (e.g. Jones et al., 1985, OG 
4.1-8). 

2) Energy independent scaling factors have been used. 

3) The resolution of the instrument is a constant fraction of the signal. 

4) Non-Z 2 corrections to the scaling laws for Zc and Zc/Zj have been determined 
directly from the data. 

Those elements contaminating the iron distribution have been approximately 
removed by scaling the contaminated distribution and subtracting according to an 
assumed set of abundances. The resulting "clean" iron distribution was smoothed and 
then scaled to the high Z elements, using cubic interpolation techniques. 

The likelihood of a given set of abundances may be calculated using Poisson statis- 
tics, and maximized by iterating until all the first derivatives are zero, yielding the best fit. 
Results from the application of this method will be reported. 
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